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1 INTRODUCTION 



ABSTRACT 

A growing number of recent observations have revealed that the Galactic globular 
cluster (GC) lo Gen is not the only GG that shows abundance spread in heavy elements 
(e.g., Fe). In order to understand the origin of the Galactic GGs with heavy element 
abundance spread ("HEAS"), we investigate the formation processes of massive GGs 
(MGGs) with masses larger than lO^M© in gas-rich dwarf galaxies interacting and 
merging with the very young Galaxy. We find that massive and compact stellar clumps 
with masses larger than 1O^M0, which can be regarded as progenitors of MGGs, can 
form from massive gas clumps that are developed through merging of gaseous regions 
initially at different radii and thus with different metallicities. Therefore it is inevitable 
that MGGs formed in dwarfs have HEAS. The abundance spread in each individual 
MGG depends on the radial metallicity gradient of the host dwarf such that it can be 
larger for the steeper metallicity gradient. For example, MGGs formed in a dwarf with 
a central metallicity of [Fe/H]= —1.1 and the radial gradient of ~ —0.2 dex kpc~^ 
can have the abundance spread of A[Fe/H] 0.2. The simulated MGGs appear to be 
significantly flattened owing to their dissipative formation from gas disks of their host 
dwarfs. Based on these results, we discuss possibly diverse formation mechanisms for 
the Galactic GGs such as M22, M54, NGG 2419, uj Gen, and Terzan 5. 
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stellar content - stars:formation 



One of remarkable recent developments in observational 
studies of the Galactic GGs is that most of the investigated 
GGs show varying degrees of chemical abundance spread: 
helium abundance spread in lu Gen and NGC 2808 (e.g., 
Bedin et al. 2004; Piotto et al. 2007), a larger dispersion in 
s-process elements abundances in NGC 1851 (e.g., Yong & 
Grundahl 2008; Milone et al. 2011), abundance spread in 
light element in "normal" GCs (e.g., Carretta et al. 2010a), 
and HEAS in M22, Terzan 5, and NGC 2419 (e.g.. Da Costa 
et al. 2009; Ferraro et al. 2009; Marino et al. 2009; Cohen 
et al. 2010). Lee et al. (2009) investigated color-magnitude 
diagrams of stars in /ifc-bands for the Galactic GCs and sug- 
gested that a significant fraction of the Galactic GGs have 
two different populations with HEAS (but see also Marino 
et al. 2009; Carretta et al. 2010b for similar and different 
claims). 

The origin of the HEAS in cj Cen has been discussed by 
theoretical models in the context of formation and evolution 
of stellar galactic nuclei in dwarfs ( "the stripped nucleus sce- 
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nario"; e.g., Bekki & Freeman 2003; Romano et al. 2010). 
Since stellar populations of stellar galactic nuclei can have 
different metallicities and ages owing to their possible for- 
mation by merging of different stellar and gaseous clumps 
(e.g., Bekki 2007), it would not be so surprising that lj Gen 
is observed to have HEAS owing to its possible origin from 
the stellar nucleus of its host dwarf. Also, it would be pos- 
sible that the formation process of other GGs with HEAS 
(e.g., Terzan 5) can be different from that of w Cen (e.g., GC 
mergers; Briins & Kroupa 2011; Bekki & Yong 2011, BYll). 
Ideally speaking, chemodynamical simulations of GG forma- 
tion based on a theoretical model are essential for discussing 
both the formation processes of GGs from gas and the re- 
sultant internal chemical abundance spread within GGs in 
a fully self-consistent manner. However no chemodynamical 
simulations on the origin of GCs with HEAS (other than lj 
Cen) have been done so far. 

The purpose of this Letter is to clearly show whether 
and how GGs with HEAS can be formed by using chemo- 
dynamical simulations of GC formation processes based on 
a GC formation scenario. We here consider that most of the 
Galactic GCs were formed in galactic building blocks (e.g., 
Searle & Zinn 1978), in particular, when they were inter- 
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acting and merging with the very young Galaxy (> 10 Gyr 
ago). We focus exclusively on [Fe/H] spread (A[Fe/H]) in 
MGCs with masses as large as or larger than IO^Mq formed 
within their host dwarfs (i.e., galactic building blocks). 



2 THE MODEL 

We numerically investigate star formation and chemical evo- 
lution processes of a gas-rich dwarf interacting and merging 
with the Galaxy in order to investigate whether MGCs com- 
posed of many stellar particles can be formed. In order to 
simulate chemodynamical evolution of the gas-rich dwarf in 
the young Galaxy, we use both the latest version of GRAPE 
(GRavity PipE, GRAPE-DR), which is the special-purpose 
computer for gravitational dynamics (Sugimoto et al. 1990), 
and high-end PCs with GPU cards (GTX 580) and CUBA 
G5/G6 software package being installed for calculations of 
gravitational dynamics. We adopt our original GRAPE-SPH 
code (Bekki 2009) which combines the method of smoothed 
particle hydrodynamics (SPH) with GRAPE for calculations 
of three-dimensional self-gravitating fluids in astrophysics. 

The dwarf galaxy is assumed to be strongly influenced 
by the fixed gravitational potential of the Galaxy in the pre- 
set study. We adopt a model of the Galaxy similar to that 
used in BYll in which the Galaxy is assumed have three 
components: a dark matter halo, a disk, and a bulge. We 
consider that (i) dwarfs forming and hosting GCs interact 
or merge with the Galaxy in the very early history of the 
Galaxy formation (> 10 Gyr ago) and thus (ii) the Galaxy 
has not yet fully developed its halo, disk, and bulge compo- 
nents at the epoch of dwarf merging. 

We assume the following logarithmic dark matter halo 
potential for the Galaxy, 

^h.io^vLiJn{r^ +<f), (1) 

where d = 12 kpc, i^haio = 131.5 km s~^ and r is the distance 
from the center of the Galaxy. The gravitational potential of 
the Galactic disk is represented by a Miyamoto-Nagai (1975) 
potential; 

GMdiak ,„N 

■l>disk = ; (2) 

y R"" + {a + VI^TV^f 

where Mdisk = 1.0 x 10"' Mq, and a = 6.5 kpc, b = 0.26 
kpc, and ii — yj -\- . The disk mass is much smaller 
than the one adopted in BYll, because we consider that 
only the first generation of the thin disk is being built when 
the dwarf starts interacting/merging with the Galaxy (> 10 
Gyr ago). We do not include the bulge in the young Galaxy, 
because we consider that the bulge can be formed much later 
from bar instability of disk stars. 

The center of the Galaxy is always set to be (x,y,z) — 
(0,0,0) whereas the initial location and velocity of the dwarf 
are free parameters that can control the orbital evolution of 
the dwarf. The initial distance of the dwarf from the Galac- 
tic center and the velocity are represented by Ri and fvVc, 
respectively, where Vc is the circular velocity at _Ri. The in- 
clination angle between the initial orbital plane and the x-z 
plane (=Galactic disk plane) is denoted as 0. The initial 
spin of the dwarf disk is specified by two angles, 9d and (pd, 
where is the angle between the z-axis and the vector of the 



Table 1. Description of the model parameters for the representa- 
tive models. 



Model no. ^ 


A4,dwil 


p c 

-^Sjdw 


kl 




J/mct X 


Ml " 


1.0 


1.8 


0.5 


-0.2 




M2 


1.0 


1.8 


0.5 


-0.1 




M3 


1.0 


1.8 


0.5 


-0.05 




M4 


1.0 


1.8 


0.1 


-0.2 




M5 


1.0 


1.8 


0.2 


-0.2 




M6 


1.0 


2.8 


0.5 


-0.2 




M7 


1.0 


1.8 


0.5 


-0.2 


0.0005 


M8 


1.0 


1.8 


0.5 


-0.2 


0.002 


HI 


1.0 


1.8 


0.5 


-0.2 





" Only the model HI is the higher resolution model with N = 10®. 
The total mass of a stellar disk in a dwarf in units of 10* Mq. 
The initial size of a stellar disk in a dwarf in units of kpc. 

The initial gas mass fraction (Afg,dw/Afs,dw)- 
^ The radial mctallicity gradient of the gas disk in a dwarf (dex 
kpc-i). 

f The value of a chemical yield (j/mct) is not shown for models 
without chemical evolution. 

angular momentum of a disk and (f} is the azimuthal angle 
measured from a;-axis to the projection of the angular mo- 
mentum vector of a disk onto the x — y plane. In the present 
study, we show the results of the models with Ri = 35 kpc, 
/v = 0.5, e = 30°, 9d = 15°, and <^d = 45°. 

We adopt the same dwarf galaxy model as that in BYll: 
the dwarf has a "cored" dark matter halo (Salucci & Burkert 
2000) and an exponential stellar disk. We describe the re- 
sult of the models with the total stellar disk masses (Ms,dw) 
being 10^ M in the present study: other models with differ- 
ent Ms,dw will be discussed in our forthcoming papers. The 
mass-ratio of dark matter halo to stellar disk in the dwarf is 
set to be 9 and thus the core radius (adm) of the dark matter 
halo is 1.5 kpc. The stellar disk has the total mass of Ms,dw 
and the size of Tis.dw and the radial (R) and vertical (Z) 
density profiles of the stellar disk are assumed to be propor- 
tional to exp(— i?/i?o) with scale length Ro = 0.2i?s,dw and 
to sech^ {Z / Zo) with vertical scale height Zo = 0.04i?s,dw , 
respectively. In addition to the rotational velocity caused by 
the gravitational field of disk and dark halo components, the 
initial radial and azimuthal velocity dispersions are assigned 
to the disc component according to the epicyclic theory with 
Toomre's parameter Q = 1.5. The vertical velocity disper- 
sion at a given radius is set to be 0.5 times as large as the 
radial velocity dispersion at that point. 

In the present study, the dwarf also has an extended 
gas disk with the total mass of Mg,dw and the size of J?g,dw 
(= 2.5i?s,dw). The disk is assumed to have an exponential 
profile with scale length 0.4iig^dw and vertical scale height 
0.04i?g,dw. The mass-ratio of gas disk to stellar one (/g) is 
an important free parameter that can control the formation 
processes of MGCs in the present study. We consider that 
the dwarf disk initially has a plenty of cold gas and accord- 
ingly an isothermal equation of state is adopted for the gas 
with a temperature (Tg) of 160 K. Star formation is modeled 
by converting the coUisional gas particles into coUisionless 
new stellar particles according to the algorithm of star for- 
mation described below. We adopt the Schmidt law with 
exponent 7 = 1.5 (1.0 < 7 < 2.0, Kennicutt 1998) as the 
controlling parameter of the rate of star formation. 
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Figure 1. The time evolution of the mass distribution of the dwarf with respect to the Galactic center (left four) and to the dwarf's 
center (right four) projected onto the x-y plane for the standard model (Ml). Gas and new stars are shown by cyan and magenta, 
respectively, and T, which is the time that has elapsed since the simulation started, is shown in the upper left corner in units of Gyr for 
each panel. The green lines in the left four panels represent the orbit of the dwarf. The locations of MGCl, 2, and 3 are indicated by 
blue, red, and green circles at T = 1.7 Gyr, respectively. 



We allocate metallicity to each gas particle according 
to its initial position: at r = _R, where r (R) is the projected 
distance (in units of kpc) from the center of the disk, and 
thus the metallicity of the gas is given as: 

[Fe/H]^^j, = [Fe/H],,^^„ + adXR. (3) 

We consider that (i) the slope ad is a free parameter and 
(ii) [Fe/H]j is determined by the total stellar mass of 
the dwarf. Accordingly we adopt [Fe/H]= —1.1 as a reason- 
able value for the adopted Af^^dw, and the observed value 
of Od ~ —0.2 for dwarf spirals (Hidalgo-Gamez et al. 2010) 
should be a reference value for the radial metallicity gra- 
dient. The chemical yield y^net and the return parameter 
Rmet are key parameters that control chemical evolution of 
dwarfs. We investigate models with (i.e., j/met > and a 
fixed i?mct = 0.3) and without (i/met = 0) chemical enrich- 
ment by supernovae. 

In order to find massive stellar clumps with masses 
(mns) as large as or larger than IO^Mq, we investigate the 
total number of new stellar particles (A'^noi) within 200pc 
around each new stellar particle. If A'^noi around a new stel- 
lar particle is larger than a threshold value nthrcs, then we 
regard the particle as being within a massive and compact 
stellar clump. Since the particle is not necessarily at the 
center of the clump, we newly determine the center of the 
clump using the position data of the neighboring particles 
around the particle. For the newly determined center of the 
clump, we estimate the total mass of the clump within the 
central 200pc. We consider that nthros ~ 1000 is reasonable, 
firstly because only massive clumps (mns ^ 5 x IO^Mq for 
/g = 0.5) can be selected and secondly because the selected 
clumps also show high stellar densities. 

The selected massive clumps are significantly larger and 



more massive than the present GCs and thus should be 
regarded as progenitors of GCs (i.e., not literally genuine 
GCs). Recent numerical studies of multiple stellar popula- 
tions in GCs have suggested that the original GC should be 
much larger and more massive than the present ones (e.g., 
D'Ercole et al. 2008; Bekki 2011). Therefore, the adopted 
assumption that massive stellar clumps can be progenitors 
of GCs is reasonable and realistic. The final massive stel- 
lar clumps are referred to as MGCs for convenience in the 
present study. Owing to the introduction of a gravitational 
softening length for new stars, the present simulations are 
not suitable for investigating the long-term internal dynam- 
ical evolution of GCs. We will discuss the long-term dy- 
namical evolution processes of GCs after their formation in 
gas-rich dwarfs in our forthcoming paper by using N-body 
codes suitable for GC internal dynamics. Table 1 summa- 
rizes nine representative models for which the results are 
discussed in the present paper. We mainly describe the re- 
sults of the "standard" model (Ml) in which /g — 0.5, 
Qd ~ —0.2, and chemical enrichment is not included (i.e., 
2/mot=0). This is mainly because we need to more clearly 
demonstrate the origin of HEAS in GCs. The total particle 
numbers for the dark matter (A'^dm), the stellar disk (Ns), 
and the gas one (A'g) used in most of the present models 
are 200,000, 200,000, 100,000, respectively. The gravitational 
softening lengths (e) for the dark matter particles, the stel- 
lar ones, and the gaseous ones are 128pc, 20pc, and 20pc, 
respectively. In order to confirm that the present results do 
not depend on the adopted numerical resolution, we also 
run a high-resolution model (HI) in which A^dm=400,000 
iVs=400,000 and iVg=200,000. 
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Figure 2. The orbits of tiie dwarf (left tliree) and the mass 
distributions of MGC2 (right three) projected onto the x-y (top), 
x-z (middle), and y-z planes (bottom) in the standard model 
(Ml). For comparison, the locations of the dwarf and MGC2 at 
T = 1.7 Gyr are shown by green and magenta circles, respectively. 
Only new stars formed from gas are plotted in the right three 
panels. 



MGC2 



Figure 3. The original location of new stars of MGCl (blue), 
2 (red), and 3 (green) within the initial dwarf disk (T = Gyr) 
projected onto the x-y plane (left) and the initial locations of 
original gas particles for new stars of MGC2 on the i?-[Fe/H] plane 
(right). The green triangle represents the final location of MGC3 
within the host at T = 3.4 Gyr (i.e., evolution into nucleus). The 
vertical dotted line in the right panel represents the size of the 
stellar disk of the host dwarf (iJs^dw)- 



3 RESULTS 

Fig. 1 shows the orbital evolution of the dwarf and the time 
evolution of the mass distributions of gas and new stars in 
the standard model (Ml). During the first pericenter pas- 
sage {Rp ~ 14 kpc) at T ~ 0.3 Gyr, the gas disk is strongly 
compressed by the tidal field of the Galaxy to form high- 
density gaseous regions. As a result of this, a number of 
massive clumps composed mostly of gas and new stars can 
be formed within the central ~ 1 kpc of the dwarf and the 
star formation rate becomes higher (~ O.OSM© yr~^). After 
the pericenter passage, a significant fraction of gas and new 
stars are tidally stripped and the stripped matter can form 
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Figure 4. The [Fe/H] distributions of new stars in MGC2 for the 
models M3 with ad = -0.05 (blue), M2 with Od = -0.1 (red), 
and Ml (i.e., the standard model) with = —0.2 (green). 



leading and trailing streams along the orbit of the dwarf 
(r = 0.7 Gyr). Some of the more massive clumps can sink 
into the central region of the dwarf owing to dynamical fric- 
tion so that they can become stellar galactic nucleus. At the 
final time step (T — 1.7 Gyr), only one massive and com- 
pact clump (identified as MGC3) with rrins = 1.2 x lO'^'M© 
(originally 1.6 x lO^'M© at T = 1.0 Gyr) can be clearly seen 
within the central 1 kpc of the dwarf. We confirm that this 
MGC3 can sink into the dwarf's center within the next ~ 1 
Gyr to become nucleus: it can be identified as an isolated 
GC only after the complete destruction of its host by the 
(older) Galaxy with the much larger mass. 

Two stellar clumps can be developed from gaseous 
tidal tails (or outer spiral arms) during the tidal inter- 
action of the dwarf with the Galaxy. Fig. 1 shows that 
the two clumps (i.e., MGC progenitors) can be stripped 
from their host dwarf and consequently located within the 
gaseous tails (T — 1.7 Gyr). These two clumps are identi- 
fied as MGCl with m„s = 7.5 x W^Mq and MGC2 with 
TTins = 1.3 X W^Mq. They can be finally regarded as "iso- 
lated" GCs owing to their large distances (> 40 kpc) from 
the dwarf. Fig. 2 shows that owing to its formation from 
the outer gas disk through gaseous dissipation, MGC2 has 
a significantly flattened shape in the x-y projection. This 
flattened shape can be seen in most of the simulated MGCs 
in the present study. 

The formation of MGCs can be severely suppressed in 
models (M4 and 5) with lower /g (< 0.2), which implies 
that MGC formation is possible only when their host dwarfs 
have a plenty of cold gas. Therefore the formation processes 
of MGCs in gas-rich disks and tidal tails described above 
can be essentially the same as those of massive stellar and 
gaseous clumps discussed in previous papers (e.g., Noguchi 
1999; e.g., Bournaud et al. 2007) in which /g is a funda- 
mentally important parameter for dynamical evolution of 
gas-rich galactic disks. MGCs can be formed even in the low 
surface brightness model M6 with larger J?s,dw (=2.8 kpc), 
which implies that mass densities of stellar disks in dwarfs 
are not so important as /g for MGC formation. 

Fig. 3 clearly shows that new stars of MGC2 originate 
from gaseous regions at different locations within the dwarf 
and thus with different [Fe/H] owing to the adopted metal- 
licity gradient of the gas disk. This result demonstrates that 
some Galactic GCs have HEAS because they formed from 
massive gaseous clumps developed from merging of different 
gaseous regions with different initial [Fe/H]. However, Fig. 4 
shows that the degree of HEAS (A [Fe/H]) depends strongly 
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Figure 5. The age-metallicity relation for MGC2 in the model 
M7 with chemical evolution (j/mct = 0.0005). Here "ages" (<agc) 
of stars are estimated at T = 1.7 Gyr (i.e., the final time step 
of the simulation): tage = means star formation at T = 1.7 
Gyr. The mean value and dispersion of [Fe/H] at each age bin are 
shown by a filled blue circle and an error bar, respectively. The 
red square represents the mean age and metallicity of MGC2. 



on ad such that A[Fe/H] is larger (~ 0.2 dex) for the steeper 
metallicity gradient (i.e., the larger absolute magnitude of 

Od). 

Fig. 5 shows that a younger population can have 
[Fe/H]~ —1 in MGC2 in the model with chemical evolu- 
tion (M7), though the mass fraction the population with 
the age less than 0.2 Gyr is only ~ 0.05. Although the pro- 
longed star formation can be seen, the vast majority of stars 
are formed in a burst at the epoch of massive clump forma- 
tion. Furthermore, there appears to be an age-metallicity 
relation that younger stars are more metal-rich in MGC2. 
It should be stressed that the relation can depend on how 
the supernova ejecta can influence the later star formation 
within cluster- forming massive gas clumps, which needs to 
be investigated by more sophisticated models in our fu- 
ture works: The [Fe/H] spread could be overestimated owing 
to the adopted instantaneous recycling approximation. The 
mean metallicities and internal [Fe/H] spread of MGCs are 
larger in models with larger Umct (e.g., M8) in the present 
study. The physical properties of MGCs (e.g., masses and 
locations) are not different between models with different 
numerical resolutions (Ml and HI). 



4 DISCUSSION AND CONCLUSIONS 

The present chemodynamical study has first demonstrated 
that MGCs with HEAS can be formed from massive gas 
clumps developed from merging of different gaseous regions 
with different metallicities. Therefore, MGCs with HEAS 
can be formed even without chemical enrichment by super- 
novae during their formation. The present study suggests 
that more massive GCs are more likely to have a larger de- 
gree of HEAS owing to their formation from gas clouds de- 
veloped from merging between a larger number of different 
gaseous regions with different metallicities. The low-mass 
GCs can be formed from single gas clouds so that they are 
likely to have no/little HEAS. The present study thus sug- 
gests that the origin of HEAS in some Galactic GCs can be 
closely associated with HEAS of ISM in their host dwarfs. 

The large metallicity spread of various elements and 
possible age variation in oj Cen (e.g., SoUima et al. 2005) 
can be consistent with the stripped nucleus scenario. How- 
ever, this does not necessarily mean that all of the Galactic 
GCs with HEAS (e.g., NGC 2419, M22, and Terzan 5) were 



formed from stripped nuclei of dwarfs. A GC with distinct 
two peaks in the [Fe/H] distributions of the stars might well 
be consistent with GC merging with different [Fe/H] (e.g., 
Briins & Kroupa 2011; BYll). A GC with a smaller degree 
of HEAS yet no clear bimodal [Fe/H] distribution might well 
form from massive gas clumps of the host dwarf and then be 
stripped from the dwarf without sinking into the center (i.e., 
without becoming the stellar nucleus) to finally become the 
Galactic halo GC. We lastly suggest that nucleation, GC 
merging, and merging of gas clouds with different [Fe/H] 
can be all promising mechanisms for the formation of the 
Galactic GCs with different degrees of HEAD. 
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